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Abstract. Recently, the MID-infrared Interferometric instrument (MIDI) at the VLTI has 
shown that dust tori in the two nearby Seyfert galaxies NGC 1068 and the Circinus galaxy 
are geometrically thick and can be well described by a thin, warm central disk, surrounded 
by a colder and fluffy torus component. By carrying out hydrodynamical simulations with 
the help of the TRAMP code ( Klahr et al. 1999), we follow the evolution of a young nuclear 
star cluster in terms of discrete mass-loss and energy injection from stellar processes. This 
naturally leads to a filamentary large scale torus component, where cold gas is able to flow 
radially inwards. The filaments open out into a dense and very turbulent disk structure. In a 
post-processing step, we calculate observable quantities like s pectral ene rgy distributions or 
images with the help of the 3D radiative transfer code MC3D jWolj2003 ). Good agreement 
is found in comparisons with data due to the existence of almost dust-free lines of sight 
through the large scale component and the large column densities caused by the dense disk. 

Key words. Galaxies: nuclei - Galaxies: Seyfert - ISM: dust, extinction - Radiative trans- 
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1. Introduction 

Within the current model of Active Galactic 
Nuclei (AGN), two observed classes of galax- 
ies can be explained by one intrinsically unique 
AGN model, the so-called Unified Scheme 
of Active Galactic Nuclei: a hot accretion 
disk surrounds the central black hole and 
crosses over into a geometrically thick torus. 
Nowadays, it is thought that the gap between 
these two i s filled by the s o-called Broad Line 
Region, see lGaskelll(l2008l) . The optically thick 
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torus blocks the sight towards the centre, when 
viewed edge-on and enables the viewing of the 
characteristic emission peak of the accretion 
disk and the Broad Line Region clouds when 
viewed face-on. Numerous observations back 
this model and also hint towards a clumpy na- 
ture of the gas and dust torus. 

2. Our hydrodynamical model 

As was recently shown bv lDavies et al. I (l2007h 
with the help of adaptive optics techniques, 
a large number of Seyfert galaxies possess 
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Table 1. Parameters of our standard model. 



Parameter 


Value 


Parameter 


Value 


Parameter 


Value 




6.6 • 10' Mq 




5pc 


T 

^ mi 


2.0 • 10^ K 


Mr, 


6.7 • 10** Mq 


cr. 


165 km/s 


M„ 


6.0 • IQ-'Mo/CyrMo) 




1.2 • 10" M© 


a 


0.5 


MpN 


2.2 Mq 




25 pc 


r 


5/3 


SNR 


IQ-'^SNe/CyrMo) 



Mass of the black hole (Mbh), stellar mass within the core radius of the distribution (Mrc), initial gas mass 
(M™J, cluster core radius (R^), torus radius {Rj), stellar velocity dispersion (cr,), exponent of the angular 
momentum distribution {a), adiabatic exponent (F), initial gas temperature (Tini), normalised mass injection 
rate {M^), mass of a single injection (Mpn) and the supernova rate (SNR). 



young nuclear star clusters. We assume that 
they are built-up in a short-duration starburst 
and are interested in the effects of the evolu- 
tion of their stars after the violent phase of su- 
pernova type II explosions is already over, af- 
ter approximately 40 Myrs. After this phase, 
energy input into the interstellar medium is 
dominated by the ejection of planetary nebu- 
lae and energy input is provided by supernova 
type la explosions. We model both of these 
processes by discrete input of mass and energy 
into a number of grid cells. The gas blobs re- 
sembling planetary nebulae additionally get a 
certain bulk velocity, made up of the turbulent 
velocity of the hot stellar system and the rota- 
tional velocity of the stars. For the rate of mass 
loss, we use the appro ximative calculations of 
lJungwiert et alJ (1200 ih . The supernova rate is 
difficult to determine and hence a more or less 
free parameter in our simulations, constrained 
by observations and a large parameter study. 
The third major ingredient of our simulations 
is optically thin cooling, where we employ an 
effective cooling curve calculated wit h t he help 
of the CLOUDY code (iFerlandl 19931 ) bv lPlewal 
(Il995h . The parameters of our standard model 
are summarised in Table [T] 

3. Resulting structure 

Within the first few percent of a global orbit 
of the simulation, single ejected planetary neb- 
ulae are visible, as well as shock fronts due 
to recent supernova explosions. The turbulent 
nature of these processes, together with cool- 
ing instability soon leads to the formation of 
long filaments of cold gas and voids of hot 



gas in-between, as evident from Fig. [Ha). Gas 
moves mainly radially inwards within the cold 
filaments, leading to a separation of gas ve- 
locities between different temperature phases. 
The largest infall velocities are reached by the 
cold material and the smallest by the hot gas 
between the filaments. This is caused by the 
larger support against gravity by thermal pres- 
sure in the hot phase and the additional turbu- 
lent pressure, arising from the perpetual stir- 
ring by supernova explosions and the input of 
mass with high random velocities. Therefore, 
cold gas moves radially inwards, until the an- 
gular momentum barrier is reached and a disk 
forms. Due to the infall of gas within the fila- 
ments from different directions, it shows turbu- 
lent behaviour. From Fig. [Ila^b) it is also evi- 
dent that density and temperature are comple- 
mentary, which results from the implemented 
optically thin cooling, which scales propor- 
tional to the square of the gas density. Pressure 
equilibrium is not reached, because of cooling 
instability in the filaments as well as recent su- 
pernova input, both acting on smaller timescale 
than relaxation processes (Fig.[T]:). The result- 
ing multi-phase medium reflects the shape of 
the cooling function as well as the local im- 
balance with energy and momentum input by 
supernovae and planetary nebulae. 

4. Comparison with observations 

With the help of a constant dust-to-gas ratio, 
we convert the resulting gas density distribu- 
tion (see Fig. [T^) into the corresponding dust 
density distribution. The latter is used in radia- 
tive transfer calculations with MC3D to calcu- 
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Fig. 1. State of the torus after roughly ten global orbits. Shown are from left to right the density distribution, 
the temperature distribution and the pressure distribution of a meridional slice through the 3D data cube. 
All images are displayed in logarithmic scaling. 



late observable quantities like spectral energy 
distributions or images at various wavelengths. 
Fig. |2(a)| shows the comparison of our stan- 
dard model at 30° inclination angle with two 
intermediate type Seyfert galaxies (scaled) ob- 
served with the Spitzer satellite. The largest de- 
viations at around 35/im are caused by crys- 
talline forsterite in emission, which is currently 
no t included in our dust model. As was shown 
by lShi et al.1 (f2006l) . AGN follow a linear rela- 
tion between neutral hydrogen column density 
and the lOjum silicate feature strength. We de- 
termine the latter as the difference between the 
simulated spectral energy distribution (SED) 
and the underlying (spline-fitted) continuum. 
Fig. |2(b)| shows the comparison of the observed 
relation for the whole Shi et al. (2006) sam- 
ple (black line) and the Seyfert galaxies alone 
(magenta) and the relation resulting from our 
hydrodynamical simulations described above. 
We obtain a shallow slope and scatter compa- 
rable to the observations. In contrast to this, our 
pr evious continuous torus m odels as described 
in lSchartmann et al] ( l2005l) are not able to ex- 
plain the full range of neutral hydrogen column 
densities without resulting in unphysically sil- 
icate feature strengths. Additionally, continu- 
ous models are in disagreement with the ob- 
served large scatter of the relation. Therefore, 
it can be taken as a further observational in- 



dication for dumpiness of AGN tori. The rea- 
son for the good comparison of our hydrody- 
namical models is caused by the two compo- 
nent structure, made up of a dense thin disk 
and a fluffy large scale torus. The filamentary 
structure guarantees the existence of more or 
less unobscured lines of sight towards the fun- 
nel region, where the silicate feature in emis- 
sion is produced and thereby relativises the sil- 
icate feature in absorption. The smearing out of 
the silicate feature in emission due to clumpy 
dust distribu tions is already a well known is- 
sue (s ee e. g. Nenkova et al. 2002, H onig et 
120061 or [Schartmann et al. 2008 ). The dense 
disk component enables also large hydrogen 
column densities. 

5. Summary 

A physical model for the built-up and evolu- 
tion of gas tori in Seyfert galaxies is presented. 
We follow the evolution of a young nuclear star 
cluster, during the phase, where the ejection of 
planetary nebulae is the dominant process for 
mass input. Together with supernova energy 
input and the onset of cooling instability yield 
a large scale, filamentary torus component with 
cavities of hot gas in-between. Within the fil- 
aments, gas cools and gets accreted in radial 
direction. Close to the angular momentum bar- 
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Fig. 2. a) Co mparison of our standar d model (solid line - ; - 30°) with IRS Spitzer observations 
of NGC4151 jWeedman et al.ll2005b and Mrk841 (H. W. W. Spoon, private communication), 
b) Comparison of resulting silicate feature strengths and column densities of our hydromodels 
(standard model plus supernova and mass loss rate study) given as s tars and the continuous 
TTM-models (standard model and models from a dust mass study, see I Schartmann et al. 
given as triangles. The filled stars denote our hydrodynamical standard model. Color denotes 
inclination angle (blue - 8°, green - 30°, magenta - 60°, red - 90°, black - intermediate angles 
of the continuous standard model). The solid lines are linear fits to the observational data shown 
in lShi et al]|2006l (black - all objects, magenta - Seyfert galaxy sample). 



rier, it opens out into a dense disk component 
with turbulent character due to the inflow of fil- 
aments and clumps of gas from random direc- 
tions. This two component structure enables a 
good comparison with observations. 
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